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I. INTRODUCTION 



Recently, the soft-collinear effective theory (SCET) has made great simplifications on the 
proof of factorization in B meson decays [l| and high energy hard scattering processes 
including resummation of large logarithms in certain regions of phase space, for example, 



e + e annihilation into two jets of thrust T 
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the deep inelastic scattering 



GDIS) in 



the threshold region x — > 1 jj] and Drell-Yan (DY) process in the case of z — > 1 ^. The 
reason for these facts is that SCET can be viewed as an operator realization of the pQCD 
analysis when the modes participating the interactions of interest are soft and collinear, 
just like chiral dynamics vs QCD at low energy region. This effective field theory (EFT) 
provides a simple and systematic method for factorization of hard, collinear and usoft or soft 
degrees of freedom at operator level, especially usoft modes can be decoupled from collinear 
modes in the Lagrangian at leading order by making a field redefinition, and the large double 
logarithms such as (a s log 2 ^r) n , where Q, A are two typical scales that characterize a process, 
can be resummed naturally through the running of renormalization group equation (RGE). 

However, all the above works have not discussed the transverse momentum (Qt) dis- 
tributions of high energy hard scattering processes. In this paper, we will investigate the 
resummed Qt distributions jfj, taking the Higgs-boson production via gluon fusion in small 
Qt region [l, 0] as an example, within the framework of SCET. It can be seen that in 
SCET the Qt resummation formula automatically separates the process-dependent Wilson 
coefficient and universal anomalous dimension of the effective operator in a process, which 
once has been studied by the authors of [^ within the Collins-Soper-Sterman (CSS) frame. 

The paper is organized as follows. In section II we start by reviewing the basic steps 
for factorization and resummation in SCET. In section III we apply it to derive the Qt 
distribution at small Qt region directly, which confirms the CSS formula. In section IV, 
we also discuss a similar formula for joint resummation. Section V contains our concluding 
remarks. The details of calculation are given in Appendix. 

II. PRELIMINARIES 

SCET is appropriate for the kinematic regions of collinear and usoft (soft) modes with 
momenta scaling: p c = (p + ,p~,p±) = (n • p,n ■ p,px) ~Q(A 2 ,1,A) and p us ~ Q(A 2 , A 2 , A 2 ) 
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or p s ~ Q(X, A, A), where A 1 is the scaling parameter, and the light-like vectors n = 
(1, 0, 0, 1), n — (1, 0, 0, —1) satisfy n ■ n — 2, and the perpendicular components of any four 
vector V are defined by V± — — (n ■ V)n'- l /2 — (n ■ V)n fl /2. 

In constructing SCET, one should first identify the scaling of all possible modes of initial 
and final states with soft and collinear degrees of freedom, then integrate other degrees of 
freedom, and the remaining modes must reproduce all the infrared physics of the full theory 
in the region where SCET is valid, which is ensured by the method of regions for Feynman 
integrals with massless quarks and gluons 1 The EFT describing the usoft(soft) and 

collinear modes is known as SCETmi), and to distinguish the two theories, the scaling 
parameter corresponding to SCET n is denoted by rj ~ A 2 , i.e., p c ~ Q{rf,l,rf) and p s ~ 

The elements of SCETi consist of usoft sectors {q us ,A us } and collinear sectors {£„, A n } 
moving in the n-direction, which are expanded as 



x 



£e-^0 niP (x), p=p + k, (1) 



where k ~ QX 2 resides in the space-time dependence of <fr n , p (x), i.e., d(p n ^(x) ~ {QX 2 )4> n , P {x), 
and p ~ Q{0, 1, A) is called label momentum, and the label operators V, V± are defined by 
picking out p~, p± momenta for collinear fields <j) n (x) 2 , respectively. The Wilson line for 
n-collinear fields has the form of 



WJx) 



eX P {^f U ■ Anix)) 



■ perms 

which is required to ensure collinear gauge invariance. The Lagrangian of collinear sectors, 
which is invariant under the usoft and collinear gauge transformation, at leading order 3 (LO) 
in A is 



£ C9 = ^{[^ + 9A£, iV v + gAl}} 2 

1 



+2Tr{c n [iP M , [iV* + A%, c n ]]} + -Tr{[zV„ A%\} 2 , (2) 

a 

£ cq = Uin -D + i fll--^- c i Pl]^n. (3) 



1 In the presence of masses, the regions analysis is very complicated, and we only discuss the massless case. 

2 The convention <fi n (x) = 4> n , P (x) for collinear fields will be used for convenience. 

3 We'll restrict our discussion only at this order through the paper. 
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Here the third line are the gauge fixing terms with parameter a and c n denotes collinear 
ghost field, and 

iT >» = V— +Vl + (in-d + gn- A us )—, 
in - D = in- D us + gn ■ A n , iD us = id + gA us , 
in-D c = V + gn- A n , %D\ =V±+ gA^. (4) 

The Lagrangian of soft sectors in SCET is identical to that of QCD. 

As for SCETn, it was emphasized that it can also be viewed as the EFT of SCETi, and 
is the final theory This suggests a short path to go into SCETn from SCETi, if the 
following matching and running steps are taken [ill ]: 

(1) Matching QCD onto SCETi at a scale /i 2 ~ Q 2 with p 2 c ~ Q 2 X 2 ; 

(2) Decoupling the usoft-collinear interactions with the field redefinitions, = Y^ffi and 
A n = Y^AfflY n . Here Y n (x) = P exp(ig J dsn ■ A us {ns + x)) is the usoft Wilson line of usoft 
gluons in n direction from s = to s = oo for final state particles, and P means path-ordered 
product, while for initial state particles, Y n is from s = — oo to s = and the daggers are 
reversed. This step leads to £ c (£ n , A m n ■ A us ) = C c (^\ Affl , 0); 

(3) Matching SCETi onto SCETn at a scale /i 2 ~ Q 2 X 2 with p 2 c ~ Q 2 rf. Thus, the soft and 
collinear modes are decoupled in the Lagrangian of SCETn- 

Next, we extend SCET to include the possibility of collinear fields moving in different 
light-cone directions n±, 712, n^, .... These directions defined by rij and Uj satisfy n^ ■ nj ^> A 2 
for i j. For simplicity we will only consider the case of head-on jets corresponding to 
collinear particles moving in the n and n directions. Since the effective theory only takes 
account the interactions of the modes in the local way, the Lagrangian of the effective theory 
contains no direct coupling of collinear particles moving in the two separate directions, 
however the usoft gluons can mediate between them in SCETi. Hence the Lagrangian in 
this case can be written by 

fc _ fC 1 fC 

*-{n,n} S ~r '-'ni 

and the soft parts are unchanged, so the decoupling transformations are also valid here. 

To illustrate the application of SCET and warm up, we consider the Sudakov effect of 
quark electromagnetic form factor in QCD Qj, i.e., the double logarithmic asymptotic of 
conservative current j M = in the following kinematics: 
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(a) nearly on-shell case 



Q 2 = -(pi - P2) 2 > -v\ = -v\ ~ A^ CD , 

Pi ~ <5(^ 2 ,i,^), P2 ~ Q(i,v 2 ,v), v ~ a qcd/Q; 



(b) off- shell case 



Q 2 = - (Pi - P2) 2 > -pi = ~p\ ~ QAqcd, 
Pl ~Q(A 2 ,l,A), p 2 ~Q(l,A 2 ,A), A ~ JAqcd/Q. 



Here pi, P 2 are the momenta of the initial and final quarks. In the above two cases we have 
omitted quark mass effects. 



Following the treatment of heavy-to-collinear current discussed in |11L [13( for case (a), 
we first match the full current onto the corresponding operator in SCETj. At LO in A, it 
gives jjj 

J" = [tnWn}rC q {V\P,li 2 )[WlU (5) 

O effective operator is de- 



lll | implies that the Wilson 



By the requirement of collinear and usoft gauge invariance, the 
termined uniquely, and the re-parameterization invariance(RPI) 
coefficient satisfies C q {V\P,n 2 ) = C q (V j ■ V^ 2 ). 

Obviously, the tree level matching condition for j 1 * leads to C q (Q 2 ,fi 2 ) = 1 + 0(a s ). 
We certainly can determine 0(a s ) correction by adopting dimensional regularization 4 (DR) 
to regulate UV and IR divergences to compute on-shell matrix elements on both sides of 
Eq.(|5]). for which is valid at the operator level and the matching calculation is independent 
of regularization method. With this choice, the fact IRqcd = IRscet provides us a direct 
matching calculation to read off the Wilson coefficients and anomalous dimensions of the 
operators in SCET. As all the on-shell loop integrals in SCET are scaleless and vanish, 
IRscet = — UVscet, and then IRqcd = — UVscet- Furthermore, the self energy diagrams 
in full QCD with massless quarks are also vanish, and all the wave function renormalization 
constants and the residues of the related propagators are equal to unity, and we also note 
that the conserved current in full QCD needs not to be renormalized. Thus, what we need 



4 The MS scheme, i.e., d — 4 — 2e and n 2 — > /i 2 e 7E /Air is used through this paper, where je is Euler's 
constant. 
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n -' 



FIG. 1: Graphical representation for quark current matching. 

to calculate is an one particle irreducible diagram, Fig^ in the full theory, which is given 
by 



Wbi) = (P2\f\pi) tree + (P2\f\pi)° ne ~ loop 
= u(p 2 )'f{l + V q )u(p 1 ), 



(6) 



V n 



a s C F .fi 2 e^. e r(l-e) 



An 

a s CF 
An 



r(l-2e) 



2 _ 3 
~2 ~~ 7 



7T 

y 



2 1, , w 2 



(7) 



Here C F = (N 2 - 1)/2N C for SU(N C ), and iV c = 3 for QCD. The e-poles in Eq.© are of IR 
character, whose opposition are just the UV poles in SCET. Thus the matching calculation 
at one-loop level gives 

a s C F . 



C q (Q\Q 
Z v 



1 + 



An 



7 (n) 

v 

in 



1 + 



s^F 

An 



~6~ 



ir 



(8) 



-(2 log ^ + 3) 



(9) 



Here Zy is defined as MS renormalization constant of the effective operator, and u has been 
set to Q to minimize the logarithms in the Wilson coefficient. It was pointed out |l| that the 
anomalous dimension of the effective operator is independent of its spin structure, for which 
can be factorized out from loop integrals. This means the evolution equation in SCET is 
universal, and only the Wilson coefficient is pro cess- dependent. 
From Eq.flSD, we obtain the RGE of C g (Q 2 ,/i 2 ), 



d\ogC q (Q 2 ,fi 2 
cflogO) 



7i(» 



dZ 



(i) 



dg 



(10) 



(a) (b) 
FIG. 2: Factorization of on-shell form factor (a) and off-shell form factor (b) in SCET. In (b), 
the soft Wilson lines are terminated and the external quarks are amputated, therefore, (b) can be 
taken as a sub-diagram of (a). Both diagrams are depicted under the gauge n ■ A n = n ■ A n = 0. 

Q2 

7i(/i) = A q (a s ) log — + B q (a s ) 

= -^<4*g + 8 ). (11) 

Here A' 1 ' = Cp and = — §CV 5 . With Eq. (fTU|) . we can resum the terms such as double 
logarithms from the scale ~ Q 2 down to the scale ~ Q 2 A 2 , we abbreviate this matching step 



as a chain QCD \qu — > SCETj 



Next, We decouple the usoft and collinear modes by the field redefinitions, which results 

in 

(pMnW n mwMpi) — (pMS ] wP }\n)r(mYnY n ]\n)(n\[w^^ ] }\Pi), (12) 

where T means time ordering operator. 

For the final step we integrate out all the off-shell modes of order ^JQAqcd and go into 
SCETn- We can rename the usoft fields as soft fields for the usoft degrees of freedom 
scaling as soft ones, and then lower the the off-shellness of the collinear fields that would be 
matched onto SCETn. Since the leading collinear Lagrangians in SCETi and SCETn are the 
same and (u)soft and collinear fields are decoupled at LO in (A)?7, all possible time-ordered 
products involve collinear fields agree exactly and we can simply replace ZgW? 

5 The notion A = ^2 n (a s /Tr) n A^ n \ we adopt {A, B} to distinguish the well known coefficients {A, B} in 
pQCD. 
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and W^tffl — > W^ 1 ^* 1 , where the superscript II denotes SCET n will be dropped from now 
on. 

Because p\,P2 of (a) are described by the collinear modes in SCETn, the general matching 
structure of SCETn diagram is shown in Fig.2.(a). The Wilson coefficient at this step is 
unity and anomalous dimension is the same as the first step, except that it runs from the 
scale ~ Q 2 \ 2 to the scale ~ Q 2 rf. We abbreviate this step as SCETi\q2 X 2 =^> SCETu\q2 T] 2. 
Collecting all the results above, we obtain the known Sudakov form factor S^(Aq C y),Q), 
leaving other coefficients omitted, 

'Aqcd 



5^(A QC D,Q)=exp(- f Q 7l (/^logA (13) 

\ -'Aqcd / 



For case (b), it can be taken as a sub-diagram of the on-shell case, from kinematical 
considerations, of which the external legs are amputated. Thus, step (1) is unchanged, and 
in step (2), (0|T[Y^y n ]|0) changes into 

poo poo 

- / / dsdte iQX (s+t \Q\T\Y n {0,ns)YU0,-nt)]\Q), (14) 
Jo Jo 

where 1 / {QX 2 ) is the effective contour length 12J and 

Y n (0,ns) = Pexp(ig £ d(3n- A us (np)^j, (15) 
Y n (0,-nt) ee Pexp (-igj^dPn-Ani-nP)). (16) 

Because of the jets Ji, J 2 with fluctuations —p\ = —p\ = QAq C d Q 2 if m Fig-2-(b), they 
must be integrated out in SCETn, and only (fT4^) is left over after step (3) associated with 
renaming the usoft modes in SCETi as the soft modes in SCETn, of which the running 
behavior is the same as Fir of jl^ . Finally, the Sudakov factor in the off-shell case is 

Si b) (Qv,Q) = exp(- / Q 7i(^)rflog^- / QA 72 (/i)rflog/i), (17) 

\ JQX JQv J 



'QX JQv 

7T (c£ Tj 



where 72 (fi) is the anomalous dimension of (|14|l . We conclude this section with a chain for 
the off-shell case, QCD\q2 — > SCETi\q2 X 2 — > SCETh\q2 v 2. Now we are ready to turn 
into the Qt resummation in the following. 
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III. METHOD OF Q T RESUMMATION IN SCET 

Since SCET is powerful to disentangle the soft and collinear interaction, and IR power 
counting ^] tells that the singular terms of Qj- distribution for DY-like processes in the 
limit of Qt — » originate from soft and collinear modes, which are emitted by partons 
from hadrons pi,P2, it is not unexpected that SCET can be applied to treat them and to 
derive the resummed part of full transverse momentum distribution for these semi-inclusive 
processes, while the remaining regular terras V flfl3 and the prescription of incorporating 
non-perturbative region (Q T ~ A-qcd) are neglected in this paper. 

For the sake of simplicity, the process of Higgs-boson production is taken as a demon- 
stration, but the method we used is not confined to this example. The dominant process 
for Higgs-boson production at the Large Hadron Collider (LHC) in the Standard Model are 
gluon fusion through a heavy quark loop, mainly the top quark, pi{P\) + ^2(^2) —* 99 - > 
(p(Q) + X with Pi = (0, 2p, 0), P 2 = (2p, 0, 0) and S = Pf P 2 + . It is convenient to start from 
the effective Lagrangian for one Higgs-boson and gluons coupling [l^ . 

C* gg = r(a s )<f>G%G^, (19) 

where r(a s ) = a ^ (y^GV) 1 ^ 2 + 0(a 2 s ) and Q = m^. Therefore, the operator for Higgs- 
boson production is 7i = G^G^f . Here the coupling a s suffers the QCD correction, which 
is unlike the case of electro-charge coupling. Furthermore, because the renormalization 
constant of ot. s G a G^ is unity up to 0(a s ), the renormalization constant of 7i is just Z~ 2 , 
where Z g is the renormalization constant of gauge coupling-^. 

If we set A 2 ~ Qt/Q with Q 3> Q T Aq C d, the situation is much like that of quark 
form factor- (a) discussed in last section, and the matching and running procedure can be 
followed. The operator 7i can match at LO in A onto 

= ~B£ a C g (Q 2 ,n 2 )B a nlx , (20) 

where 

K = M?, 

<?r = + gA», iV v n + gA%W n , (21) 
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P2 



Pi 





xC g (Q 2 ,fM 2 



FIG. 3: Graphical representation for gluon current matching, 
and n <-> n for Q£ v . 

The one loop calculation, FigJHl is similar to quark current, except for dividing the final 



result by Z 2 . Finally, 



[9i\n\g 2 ) = (gi\n\g 2 ) tree + ( gi \H\g 2 ) 



i one— loop 



+ C.t. 



( gi \H\g 2 ) tree (l + V g ), 



(22) 



V n 



An 
Air 



r(l-e) /^ 2 e 



2C A _ 2A) ^j/ N 



r(l -2e) \ Q 2 J \ e 2 e 
- 2 -^--(2C A \o g ^ + 2(3 ) 

2 P 2 no , /t h i Cavt 2 



-C A log^-2/3 log^ + ^ + 



is obtained 6 , from which we read 



C 9 (Q 2 ,Q 2 ) = l + ^-(Af + ^f 



An) 
J H 



An 



^ + i(2C A log^ + 2/3 ) 



(23) 

(24) 
(25) 



Here C A = N c , flo = ^Ca — |^/Tr, A g = 11 + 2n 2 , Tr — \ and nf = 5 is the number of 
active quark flavors. Then the RGE of C g (Q 2 ,/j, 2 ) is 



dlogO*) 



<9# 



(26) 



Q 2 



7i(/i) = A 3 (a s ) log — + B g (a s ) 

P 1 

= -(^)(C A log^ + /3 ), 



(27) 



with = and = —(3q. Thus the evolution from the scale ~ Q 2 to the scale 
C g (Q 2 ,Q 2 X 2 ) = C g (Q 2 ,Q 2 )exp (- £^ ^7i(a*))- (28) 



~ <5 2 A 2 gives 



Here we have absorbed the scale dependence of a s into that of Z„. 
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As shown above, the extraction of A, B in SCET is different from that of A, B in pQCD, 
i.e., there is no need to calculate real correction which is more difficult to handle. Using the 
virtual part of higher order calculation, such as the two loop on-shell quark and gluon form 
factor |l6j], we can find the 0(a 2 ) universal anomalous dimension. For example, 

Ai 2) = \c a K, K = C A (% - y) - ™T R n f , (29) 

where C q = Cp, C g = 

After performing field redefinitions, Eq. (|2fj|) can be directly used to running Ti. from the 
scale ~ Q 2 to the scale ~ Q 2 f] 2 without loss of degrees of freedom. So the relevant operator 
for Higgs-boson production at the scale ~ Q 2 V 2 {Qt) * s 



H = C g (Q 2 ,Q 2 V 2 )Tr{T[Y^Y n B^Y n B n ^}} 

1, 
2 

= C g (Q 2 ,Q 2 V 2 )H. (30) 



c 9 (Q\Q 2 v 2 hnyfy a n^ c n,} 



Here y n tn) is the adjoint soft Wilson line from — oo to in n{n) direction for incoming fields. 
Now, we have completed the procedures corresponding to step (1), (2) and (3). 

To obtain the differential cross section, we relate it to the composite operator 7i at the 
renormalization scale /i 2 ~ Q 2 rj 2 in SCETn, where the cross section can be written as 

1 da Icsum 



afj dQ 2 dydQ 2 T 

->2M2 



Hj(Q)e- s "^V S cET(QT, Q, H\ (31) 



where Hj(Q) = \C g (Q 2 ,Q 2 )\ 2 is a function of a s (Q), and 

r(0) . . u /o n y s {QWu 



- (V2G F )^^5(Q 2 -m 2 H ), (32) 
Q 2 d{i 2 r " ^ 2 



Q 

A g (a s )log— + B g (a s 



(33) 



and 



1 Ax SCET (/i) 

ascMQT ' Q ^ ) = 4$dQ 2 dydQ 2 T (34) 

represents the normalized differential cross section calculated in SCETn with the composite 
operator Ti.. The general structure of relevant diagram is shown in Fig. where the soft 
and collinear modes are decoupled and the spin and color are summed over in the matrix 
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FIG. 4: General structure of SCET cross section for Higgs-boson production under the gauge 
n ■ A n = n ■ A n = 



element of hadron, from which the SCET cross section can be written in the form of multiple 
convolution, 

o"scet(<2t, Q,y) = J d 2 k T1 d 2 k T2 d 2 k TS 5 2 (k T1 + k T2 + k TS - Q T ) 

xj pi (xt, kri, fi)J P2 (x 2 , k T2 , n)S(k T s,n), (35) 

where x x = Qey/^fS, x 2 = Qe~y/^S for Q\ < Q 2 , and 

J pl (x u k Tl ^) = -^j^-rr, I dy + d 2 y ± e-^ p -y + -^y^ 
X\F \Itx) 6 J 

x(p 1 |Tr[^(y + ,0,^)i3 nQ (0)]b 1 ), (36) 
J P2 (* 2 , W) = ^T^p / dy-d 2 y x e-^ + y-^ 

x (p 2 |Tr[£?(0, y~, y ± )B na (0)]\p 2 ), (37) 



X' 



T[J^ffl(0)|n>, (38) 



with T denoting the anti-time ordering operator. Obviously, in Eq.(38) the matrix element 
has been factorized, and the delta function is imposed by momentum conservation. 

Next, to factorize the phase space, the trick of Fourier transforming to impact parameter 
space is significant 

J d 2 Q T e^ T 5 2 C£ k Tl - Q T ) = J] e* fc «. (39) 

i i 

Then, for each transverse momentum kn, one obtains 

J d 2 k Tl e A ^f(k Ti ) = f(b). (40) 
12 



This produces the simple product 

5-scetO,Q,aO = J Pl {xt,b, fx)J P2 (x 2 ,b, fx)S(b, fx). (41) 

Because of KLN theorem, the contributions from the soft modes are free of IR divergences. 
So only the collinear divergences are survived, therefore after matching the SCET cross 
section onto a product of two parton distribution functions (PDFs) given by j^J, which are 
equivalent to the conventional PDFs f a /pi(xi, aO a ^ LO in A, the remaining IR divergences 
can be absorbed into these nonperturbative inputs, of which the evolutions are controlled 
by the DGLAP equations. This leads to J(| 

J Pi (xi,b,/i) = J2(fa/ Pl ®Cga)(xi,b,y) 

a 

^Hf^fa/ptti^a&Ari, (42) 
a Jx i S S 

a. 



c = Vf— ) n c {n) 

^ga — 2_^\a ) ^ga 

n=0 



If we define 



C ga (z, ^,/i) = c ga (z,b,y)[S(b,fi)] l 2, b = 2e-T* (43) 

then 

0"SCEt(&,Q,A*) = (fa/pi ® C ga )(xx, y, A*) (/ft/pa ® C fl6 ) (x 2 , y,A*)- (44) 

Obviously, C^(z) = 5 sa 5(l — z), and we have derived in Appendix that 



'ga W ga 



C«(*) = -2Pl(z) - C A ^-6 g J(l - z), (45) 



where fx has been set to bo/b to eliminate large constant factors in C^-J(z,fx), and Pg a (z) 
represent the 0(e) terms of the DGLAP splitting kernels. 

Combining Eq.([31 |) -(|44 |) and Fourier transforming back to Qt space, we obtain the re- 
sumed formula of transverse momentum distribution for Higgs-boson production in SCET, 

J_ ^ rCSUm = uHQ) r -b,UbQ T )Y e~ s ^ 
o 9 «}dQHydQ\ ^ [Q) J 2ir bJoWT) Y 

X(/a/pi ® C ga )(xx, j)(fb/ P2 ® C 9b )(x 2 , y). (46) 
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The similar reasoning leads to the general form for Qt resummation, 
1 da rcsum „ P/ „, db. 



JU * n ab 

x(fa/ Pl ® C ca )(xi, b -j-){fb/ P2 ® Ccb)(x 2 , y), (47) 



where F and c stand for the type of process and of parton participating the elementary 
sub-process, for example, F = DY and c = a for Drell-Yan process. The formula (}4T)) is a 
little different from the known CSS formula 



x^®^)^,^/^®^)^^). (48) 



Transforming Eq. (j47|) into the form as Eq.([48|) by the identity 



Hf(Q)=exp 



Q 2 dfi 2 .d\ogU F 



-(3(a 



c 



'bl/b 2 /i 2 d\oga s 
with (3{pt s ) denoting the QCD /3-function, one finds 
1 da rcsum r°° db 



Hf(6d/6) (49) 



x(/a/pi ® C'ca)^!, j)(fb/p 2 ® C a )(x 2 , y). (50) 

The corresponding coefficients {A, 5, (7} are related to {A, B, C} through equation 

A(a s ) = A(a s ), 

B{a s ) = B(a s ) - P(a s ) d ^ gU * , (51) 

d log a s 

C ab (z) = C a6 (2)[H£(6 /#, 

Up to next-to-leading-logarithmic-order (NLLO), {A^ 2 \ B«, C^} and {A( 2 \ flW, C^} is 
compatible with each other, and further calculation and confirmation are required at higher 
order. 

It can be seen that SCET provides a natural framework of Qt resummation by conven- 
tional RGE in EFT. We have noted that the matched effective operator is determined by 
collinear and soft gauge invariance and is unique. In addition, the corresponding anomalous 
dimension is independent of its spin structure and is universal, while the process-dependent 
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quantity resides in the Wilson coefficient. Even more, the coefficients C a b defined in SCET 
are process-independent too. So the matching and running procedure in EFT naturally 
separated the process-dependent and universal contributions to a process, i.e., {A,B,C} 
are universal and only Hf is process-dependent. In pQCD, the formula and relation like 
Eq . (|47|) and Eq.(|51|) have been proposed by the authors of j9(. 

Compared with SCET, pQCD analysis invoke gauge invariance and a new evolution equa- 
tion j(| which comes from differentiating the jet-function from the factorized cross section 
with respect to the axial parameter in axial gauge to separate soft and collinear contribu- 
tions, which is crucial to resum the double logarithms, since RGE in full theory only resums 
single logarithms between two scales. 

We conclude this section with a chain for the Qt resummation in this section, 

QCD\q2 — > SCETj\q2 X 2 SCET ii \q2 v 2 < — DGLAP\^2, 

where the last arrow indicates that the PDFs used at the scale Q 2 r] 2 can be obtained from 
those at some fixed scale /xjj by the evolutions of the DGLAP equations. 



IV. DISCUSSION 

(I) Applying the formula to the production of lepton pair via virtual photon, one 
can find, 

b b 
x(fa/ Pl ® C qa )(xx, j)(f b / P2 ® C qb )(x 2 , j), (52) 



where 



fQ 2 dfl 2 r " 2 



(0) „2 ^ n "em m DY (n\ 1 i uts ^ F t q i ' 2^ 



by* /i 2 



Q 

A 9 (a s )log— + B q (a s ] 
/i 



Cg } = S qa 8(l - z), C« = -2P* a (z) - C F ^8 q J(l - z) 



(53) 



(II) The reformulation of joint resummation can also be made straightforwardly in SCET. 
In fact, threshold resummation [^] under z = Q 2 /S — > 1 for dcr TCSnm /dQ 2 is performed in 
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moment-iV space, and the relevant A 2 = l/iV~l — z with N = e lE N. The conclusion of 

n 

|5J] can be represented by a chain, 

QCD\ Q 2 — > SCETj\ Q i X 2 SCET U \ QW <— DGLAP\ (l *. 

We observe that the two chains for Qt and threshold resummation in SCET have identical 
structure. This suggests that we can do threshold and Qt resummation for da resnm /dQ 2 dQ T 
simultaneously. The relevant A 2 ~ l/x(iV, 6) with b = bQ/bo is an interpolation of A 2 ~ 1/N 
and A 2 ~ 1/6, let us say \ul\ 

*^ 5) =» + TW '4 (54) 

which approaches to TV for 6 <C iV and to 6 for 6 ^> N, respectively. The matching steps for 
joint resummation then can be written as 



QCD\ Q 2 — ► SCET z \ Q 2 /x => SCET n \ Q 2 /x 2 < — DGLAP\ 

which leads to similar result as Eq. (j47p corresponding to that of Q], and the Mellin trans- 
formed and jointly resummed cross section follows, 



1 drr Tesum (M) r°° rib 

x (c m (N)f a/pi (N, Q/x)) (c- cb (N)f b/p2 (N, Q/ X ) ) , (55) 
where 4>(N) = Jq 1 d££ N </>(£) for any function 0(£) with < £ < 1 is used. 



, -S c (Q/x,Q) 



V. CONCLUSION 



We have presented the method of Qt resummation in the framework of SCET and 
given a simple correspondence between {A, B, C} in SCET and the well known coefficients 
{A,B,C} in pQCD, with which the available information is compatible. The equivalence 
of the two framework can be confirmed by higher order computation. We have also shown 
that the reformulation of joint resummation can be performed in SCET directly. So any 
process, which is confined to the soft and collinear regions by dynamics or kinematics, can 
be treated in SCET following the steps outlined above. 
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Appendix. A 

In this appendix, the details of the calculation to extract in SCET are given explicitly. 
Because is related to the emission of a soft or collinear gluon, and the phase space 
is already factorized in this case, we will exploit a special form at 0(a s ), for which there is 
no need to cover the non-perturbative region (Qt ~ Aqcd), 
1 da iesum d 



H*(Q)- 



a^dQ 2 dydQ 2 T cK ^'dQ\ 
ascET(QT,Q,/i) - j Q dq T ^ dQ2d d 2 



e~ s ^a SC ET(QT,Q,v) 



(56) 



<r(°) dQ 2 dydq\ 

= J P i{xi,Q T , fj)Jp 2 {x 2 ,Q T , v)S{Q T , y). (57) 
The same reasoning as in section 3 leads to 

J pl {x 1 ,Q T ,ii)J P2 (x 2 ,Q T ,fi)S(Q T ,fi) = ® C m )(xi,(3 T ,/i) 

a. b 

x(fb/p 2 ®Ccb)(x2,QT,y). (58) 
Then we get the differential form for Qt resummation, 

o^dQHydQl cW> dQ 2 T ^ 6 

x(fa/ Pl ® C ca )(x 1 ,Q T )(f b/p2 ® C 5b )(x 2 ,Q T ). (59) 

Here /i = Qt is set to minimize large factor in C^. Taking the Q\ integral of Eq. (j47j) and 
Rq-lpSjl from to Q\ and then expanding them at 0(a 8 ), we find C (1) (z, fy) = C (1) (z, Qt), 
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V 



= I 



yt n-p 
2 p 



a k v,b 



= -ig^S ab jz 



p 



2 V A* n-p 



n-q 



n-pn-q r 



= -gT a ^ 

•J n-k 



"dnsiroTsip 
k 



-gT a - 



2l 

n-k 



iQ. 



FIG. 5: Feynman rules for ra-direction collinear particles in Feynman gauge a = 1 UJ|. Here 
the direction of gluon momentum in the Wilson line is along the collinear particle. The rules for 
soft fields and collinear gluon are the same as that for QCD, and n <-> n for redirection collinear 
particles. 

thus we will use Eq. fJ57j) -(|58 jl to calculate C^. However, it should be emphasized that this 
formula is valid only with {A^, BOO, CCO}, ancx the two formulas must be equal at this 
order, which suggests a way to adjust the parameters in b* prescription Q]. Previously, the 
authors (DDT) of Q| have derived a similar formula, which is corresponding to our result 
at leading-logarithmic-order (LLO). Later, the extended DDT formula in the CSS frame 
with {A^ 2 \ B&\ C^} was suggested in 20j, whose coefficients {A, B,C} at 0{a s ) are just 
our {A^ 1 ', B^, C*- 1 - 1 }. The real radiative contribution to the differential cross section at the 
scale Q\ is 



do, 



SCET 



,47T, 



dQ 2 dydq\ 



\M^( Pl ,p 2 ,k,y)\ 
8s(2vr) 2 



5 1 (it Umin){'U' ^max) 
U \U 



(60) 



Here M-ab(pi,P2,k,y) denotes the corresponding matrix element, and the usual invariants 



18 



11 



n •' 




FIG. 6: Cut diagrams for the emission of gluon in SCETn, and the mirror diagrams are not shown 



are defined as 



s = ( Pl +p 2 ) 2 , t = (pi- k) 2 , 

Q 2 



z 



u = (p 2 - k) 2 , 

s 

and the two roots of the equation (pi + p 2 — q) 2 = are 

,z- 1 - 



u, 



u. 



Q 



l-z) 2 - 4zq 2 /Q 2 



2z 



Q 5 



.z-1 



1 - zf - Azq 2 /Q 2 



2z 



The expression |-M^(pi,£>2, k, y)\ 2 in SCET can be written down by cut diagrams, and 
the Feynman rules we need are shown in FigJSJ For example, in FigEl the real contribution 
of the collinear gluon in the n direction is 



\Mf g ^{ Pl ,p 2 ,k,y)\' 



4na s fi —Cf{ Tr 



S7m 



jhn ■ (pi — k) fi 

2 (pi - ky 2" 



(61) 



f (A- A± _ (A- j)± A± _ \ A 

\ u n ■ (pi — k) n ■ pi n ■ (j>i — k)n ■ p\ ) 2 

» 1 AKA- A± A± (A- /fe)± A± -A A 

2 V n ■ pi — k n ■ p\ n ■ [pi — k)n ■ pi J 2 

J}v_ fi_ fi_ v fa n-(pi-k) ^/i 
S ^n-k2 2 U 2 {pi-kf 1 2 
(pi - k)± A± 



n-jpi-k )_ M fi 
(Pi - k) 



2 r 



2Tr 



4vra s /i ie — Tr[ sl 'Yir}2C F \ (1 - e 



36 



n ■ (pi — k)n ■ pi 



(n ■ (pi — k)) 2 (n ■ pi) 2 
If we drop the common factor 



n ■ (Pi - k ) 1 2 2 n-{pi-k) 



(pi - k) 2 



n ■ k(pi — k) 2 



(62) 
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and use the following parametrization for momenta pi,P2, k, 

Vx = (o,p-,o), p 2 = (P + ,0,0), 

k = ( n q2 \ P ^ C 1 " Z ^ P ^ ( 63 ) 

1-Zi 

u = — (1 — z±)s, tu = sq?, (64) 

the above expression can be simplified to 

\Mir(PuP^y)\ 2 - 2 ^ 1(1 " £)(1 2 "" l)2 + 2Zl] . (65) 

zq T 

The contribution of n collinear gluon is given by n <-> n and £i 22, 

l*CW.MI a - 2 ^ [(1 " e)(1 r^ 2 + 2Z2] , (66) 
and the soft gluon contribution is 

AC 

\M s q f\p 1 ,p 2 ,Ky)\ 2 ^^. (67) 

The condition that the emitted gluon is to be collinear is (1 — Zj)P~ ^> gx or — > 2 
for z = 1,2, so only half of the phase space is covered, i.e., u = u min for u = —(1 — Zi)s, 
under which we could safely make the substitution Zi —> z. The soft gluon is guaranteed by 

(1 - z{)P~ ~ q T or z — > 1. Note that |.M ci(0) | 2 — > |.M s(0) | 2 when ^ — > ^ — >■ 1, which results 

• 7 
m 

l-M^Cpi.pa.fc.y)! 3 ' = l^fbi,^,^^)! 2 ^^ 

+I^ 0) (pi,p 2 ,a ; , 2 /)| 2 _^ 

« p- 

2CV[(l-e)(l-z) 2 + 2z] 
_ 2(1 -z)P m (z,e) 



zq\ 



(68) 



similarly, 



zq T 



7 We have divided soft contribution into two parts. 
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where 



C F 

C F 
2C A 



1 + z 2 



1-z 

l + (l-2) 2 



e(l-z) 



— ez 



z 



1 



1 — z z , . 

+ - + z(l-z) 

z 1 — z 

2z(l - z) 



(69) 



1-e 



These are the corresponding results of 

Next we match the Q F integral Eq. tJoTj) onto PDFs after making Mellin transformation, 



E& } (n) 



dzz % 



dQ 2 dydq 2 T ' 



1 



2q T 



(l-2) 2 -4^ 



(70) 
(71) 



where the upper limit 2 = 1 — 2q F /Q is to make the integrand meaningful, and the integral 
(171)1 can be evaluated easily if we retain only the singular terms as q\ — > 0, 



2q T 



dz 



2(1 -z)P qq (z,e 



+ / dz(z n -l)2P qq (z,e) 



Q 2 



2C F log \ - 3C F + 2 7w (n) + 2efjn), 
q T 



(72) 



where ■y qq (n) and 7 L(n) are the moment of regularized splitting function and of P qq (z). So 



dq%T!§(n) 



1 [Qt dq\ n 2 e^ E . 



r(l-e)io q T K q T ' 
2C F 1 



+ 7 (3C F + 2C F log^ 



2CV log ^ - 3C F + 2 lqq (n) + 2e 7 |, (n) 



6 



(73) 



Q 2 ,„ /i 2 



2 2 2 

+CV log 2 ^ + 3C F log ^ - 2 7( » log ^ - 2CV log ^ log 



The virtual correction comes from the UV renormalization constant 

,2 



25Z V = £ 



2 -^f + -{3C F + 2C F \og^) 
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which cancels the first two e-poles. The remaining is cancelled by the renormalization of 
PDF, i.e., 



f q (n) = / ff („) = l + i ± 
where \i — jip is implied. Finally, we get 



27 9 »(— -— ) 



cuv ClR 



C$(n,Q T ) = -27' (n) 



Similarly, 



C F 7T 2 

6 

C F vr 2 



- 2 W " 



(J(l-z), 



-2P£(z)-^W(l-*)- 



(74) 

(75) 
(76) 

(77) 
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